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Abstract. We discuss pycnonuclear burning of highly exotic atomic nuclei in deep crusts of neutron stars, 
at densities up to 10 13 g cm -3 . As an application, we consider pycnonuclear burning of matter accreted on 
a neutron star in a soft X-ray transient (SXT, a compact binary containing a neutron star and a low-mass 
companion). The energy released in this burning, while the matter sinks into the stellar crust under the 
weight of newly accreted material, is sufficient to warm up the star and initiate neutrino emission in its 
core. The surface thermal radiation of the star in quiescent states becomes dependent of poorly known 
equation of state (EOS) of supranuclear matter in the stellar core, which gives a method to explore this 
EOS. Four qualitatively different model EOSs are tested against observations of SXTs. They imply different 
levels of the enhancement of neutrino emission in massive neutron stars by ( f) the direct Urea process in 
nucleon/hyperon matter; (2) pion condensates; (3) kaon condensates; (4) Cooper pairing of neutrons in 
nucleon matter with the forbidden direct Urea process. A low level of the thermal quiescent emission of 
two SXTs, SAX J1808. 4-3658 and Cen X-4, contradicts model (4). Observations of SXTs test the same 
physics of dense matter as observations of thermal radiation from cooling isolated neutron stars, but the 
data on SXTs are currently more conclusive. 

PACS. 97.60. Jd Neutron stars - 26.60.+C Nuclear matter aspects of neutron stars 



1 Introduction 

Neutron stars (NSs) are the most compact stars known in 
' the Universe. Their masses are M ~ 1.4 Mq (Mq being 
, the solar mass) while their radii are R ~ 10 km. They 
are thought to consist (e.g., Ref. 0E|) of a massive dense 
' core surrounded by a thin crust (of mass < 0.01 M and 
thickness < 1 km) . The crust-core interface occurs at the 
density p « po/2, where po — 2.8 x 10 14 g cm" 3 is the 
standard density of saturated nuclear matter. 

The NS crust is composed of neutron-rich atomic nu- 
clei, strongly degenerate electrons and (at p > 4 x 10 11 
g cm -3 ) of free neutrons dripped from nuclei. NS cores 
contain degenerate matter of supranuclear density (up to 
10 — 15 po)- Its composition and equation of state (EOS) 
are model dependent (being determined by still poorly 
known strong interactions in dense matter). On the other 
hand, these properties are almost not constrained by ob- 
servations (e.g., Ref. U). The NS core is usually subdi- 
vided into the outer core (p < 2 po) and the inner core 
(p > 2po). The outer core is available in all NSs, while 
the inner core is present only in massive, more compact 
NSs. The outer core is thought to be composed of neu- 
trons, with the admixture of protons, electrons (and pos- 
sibly muons). The composition and EOS of the inner core 
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is still a mystery. It may be the same composition as in 
the outer core, with a possible addition of hyperons. Al- 
ternatively, the inner core may contain pion- or kaon con- 
densates, or quark matter, or the mixture of these compo- 
nents. Another complication is introduced by superfluidity 
of neutrons, protons and other baryons in dense matter. 
Superfiuid gaps are very model dependent (e.g., Ref. |3]). 

It is unlikely that the fundamental problem of the EOS 
of supranuclear matter can be solved on purely theoretical 
grounds. However the solution can be obtained by compar- 
ing theoretical models of NSs with observations. We will 
mention two lines of such studies. The first one is based on 
theory and observations of cooling isolated NSs (e.g., Refs. 
4.5 ) ■ These studies, carried out over several decades, 
have constrained the properties of dense matter but have 
not solved the problem (Sect. 4). Here, we focus on an 
alternative method to test the same physics on another 
class of objects, transiently accreting NSs in soft X-ray 
transients (SXTs). This method is based on pycnonuclear 
burning of accreted matter. The method is new but seems 
to be currently more restrictive. 



2 Deep crustal heating 

Nuclear fusion reactions in normal stars proceed in the 
thermonuclear regime, in which the Coulomb barrier is 
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penetrated owing to the thermal energy of colliding nu- 
clei. Here, we discuss another, pycnonuclear regime where 
the Coulomb barrier is penetrated due to zero-point vi- 
brations of nuclei arranged, for instance, in a lattice. The 
thermonuclear regime is realized in a rather low-density 
and warm plasma, whereas the pycnonuclear regime op- 
erates at high densities and not too high temperatures. 
Pycnonuclear reactions are almost temperature indepen- 
dent and occur even at T = 0. They were suggested by 
Gamow 8 in 1938. The first calculations of pycnonuclear 
reaction rates were done by Wildhack [§] upon Gamow's 
request. The strict approach for calculating these rates 
was formulated by Salpeter and van Horn JOj who studied 
also three other regimes of nuclear burning in dense mat- 
ter (the thermonuclear regime with strong enhancement 
due to plasma screening effects; the intermediate thermo- 
pycno nuclear regime; and the thermally enhanced pyc- 
nonuclear regime). The pycnonuclear regime has been an- 
alyzed later in a number of publications (e.g., |1 11112) and 
references therein) . Pycnonuclear reactions are extremely 
slow at densities p typical for normal stars but intensify 
with increasing p. For example, carbon burns rapidly into 
heavier elements at p > 10 10 g cm" 3 . 

Pycnonuclear burning can be important in transiently 
accreting NSs. The heat released due to the infall of mat- 
ter and thermonuclear reactions in the surface layers is 
radiated away by photons from the NS surface and can- 
not warm up the NS interiors. The accreted matter sinks 
into the NS crust under the weight of newly accreted mate- 
rial. The density gradually increases in a matter element, 
and the nuclei undergo transformations - beta captures, 
absorption and emission of neutrons, and pycnonuclear re- 
actions. The nuclei evolve then into highly exotic atomic 
nuclei which are unstable in laboratory conditions but sta- 
ble in dense matter. The transformations and associated 
energy release have been studied by Haensel and Zdunik 
[131114) . At p > 10 9 g cm -3 the transformations are almost 
temperature independent but depend on the composition 
of matter at p ~ 10 9 g cm -3 . The main energy release 
occurs at densities from about 10 12 to 10 13 g cm -3 , sev- 
eral hundred meters under the NS surface, in pycnonuclear 
reactions. For iron matter at p ~ 10 9 g cm -3 , these reac- 
tions are Q3| 34 Ne+ 34 Ne^ 68 Ca; 36 Ne+ 36 Ne^ 72 Ca; and 
48 Mg+ 48 Mg^ 96 Cr. The total energy release is then wl.45 
MeV per accreted baryon; the total NS heating power is 
determined by the mass accretion rate M: 

L dh = 1.45 MeVAf/rTiN 

w 8.74 x 1O 33 M/(1O" 1O M yr" 1 ) ergs" 1 , (1) 

where ton is the nucleon mass. It produces deep crustal 
heating of the star. This heat is spread by thermal con- 
ductivity over the entire NS and warms it up. 



3 Thermal states of soft X-ray transients 

It is possible that the deep crustal heating manifests itself 
in NSs which enter SXTs, compact binaries with low-mass 



companions |15| . These objects undergo the periods of out- 
burst activity (days-months, sometimes years) superim- 
posed with the periods of quiescence (months-decades). 
An active period is thought to be associated with an ac- 
cretion of matter from a companion through an accretion 
disk. The accretion energy released at the NS surface is 
so large that a SXT looks like a bright X-ray source with 
the luminosity Lx ~ 10 36 — 10 38 erg s _1 . The accretion is 
switched off or suppressed in quiescent periods when Lx 
drops below 10 34 erg s" 1 . In many cases, the spectra of 
quiescent emission contain the thermal component, well 
described by NS hydrogen atmosphere models with effec- 
tive surface temperatures from few 10 5 K to ~ 10 6 K. This 
can be the radiation emergent from warm NSs. 

The suggestion to interpret the quiescent radiation 
spectra of SXTs with hydrogen atmosphere models was 
put forward by Brown et al. fTSI who assumed also that 
NSs in SXTs could be warmed up by the deep crustal 
heating of accreted matter. It is important that the crustal 
heat is partly radiated away by neutrinos from the entire 
NS body. Because the neutrino luminosity depends on the 
NS structure, the remaining heat, diffused to the surface 
and radiated away by photons, becomes also dependent 
of the NS structure. This opens an attractive possibility 
(see j^.19,20,21] and references therein) to explore the 
NS structure by comparing the observed quiescent ther- 
mal radiation from SXTs with theoretical predictions. 

We outline the results of such studies using the theory 
of thermal states of transiently accreting NSs described in 
Rcf . [22] • These stars are thermally inertial, with thermal 
relaxation times ~ 10 4 yr |19| . In the first approximation, 
a transiently accreting NS is in a (quasi) stationary steady 
state determined by the mass accretion rate M = (M) 
averaged over thermal relaxation time-scales. Typically, 
(M) ranges from 10" 14 to 10~ 9 Mq yr -1 and does not 
increase noticeably the NS mass during long periods of 
NS evolution. Quasistationary states are expected to be 
rather insensitive to variations of the accretion rate asso- 
ciated with a sequence of active and quiescent states. 

Thermal states of accreting NSs can be greatly affected 
by the neutrino emission from NS interiors (mainly from 
the cores). The neutrino luminosity of the star, L u , in its 
turn, can strongly depend on the NS mass (as schemati- 
cally shown in Fig.flJ. Low-mass NSs possess nucleon cores 
with not too high neutrino emission determined mainly by 
the modified Urea (Murca) and nucleon-nucleon brems- 
strahlung (brems) processes. High-mass NSs, in addition 
to the outer nucleon cores, possess the inner cores whose 
composition and EOS are unknown (Sect. 1). Their neu- 
trino emission can be enhanced with respect to the emis- 
sion of low-mass NSs. Four qualitatively different enhance- 
ment mechanisms are (e.g., Ref. |7|): (1) a very powerful 
direct Urea (Durca) process in the cores of massive NSs 
containing nucleons (and possibly hyperons); (2) a less 
powerful direct-Urca-like process in pion-condensed cores; 
(<?) even less powerful similar processes in kaon-condensed 
or quark NS cores; (4) even weaker enhancement of the 
neutrino emission owing to mild Cooper pairing of neu- 
trons in the nucleon inner NS cores with forbidden direct 
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Fig. 1. A sketch of the NS neutrino luminosity L v versus NS 
mass at the internal stellar temperature T = 3 x 10 8 K for four 
scenarios of NS structure (from Ref. 0). 



Urea process A transition from the slow neutrino 

emission in low-mass NSs to a fast emission in massive NSs 
with increasing M occurs in medium-mass NSs; the mass 
range of medium-mass stars is model-dependent. The NS 
neutrino luminosity is a strong function of the internal 
stellar temperature. 

A thermal state of the star is determined by the ther- 
mal balance equation (which implies that the crustal heat- 
ing power is equal to the sum of photon and neutrino lu- 
minosities). Solving this equation, one obtains a heating 
curve, the dependence of the thermal photon luminosity 
(as detected by a distant observer) on the mean mass ac- 
cretion rate, L?°(M). 

Several examples are presented in Fig. 2. The upper 
curve is calculated for a low-mass NS (M — 1.1 Mq) 
whose core is composed of neutrons, protons and elec- 
trons. In the core we use a moderately stiff EOS proposed 
in Ref. JTj (the same version as employed in Ref. 0]). It 
is assumed that the core contains strongly superfluid pro- 
tons. Superfluidity suppresses the modified Urea process, 
the star has a low neutrino luminosity due to neutron- 
neutron bremsstrahlung and stays relatively warm. The 
lowest curve in Fig. 2 corresponds to the maximum-mass 
NS (ill = 1.977 Mq) with the same EOS. This star has a 
massive inner core where the direct Urea process operates 
and proton superfluidity dies out. The neutrino luminosity 
becomes exceptionally high, and the star is very cold. The 
next two curves above the lowest one are schematic mod- 
els [22] of high-mass NSs with pion-condensed or kaon- 
condensed cores. The next curve is calculated for a high- 
mass (M = 2.05 Mq) star, where the direct Urea process is 
forbidden but the neutrino emission is enhanced by mild 
Cooper pairing of neutrons in the inner core (the same 
EOS [23 an d model for superfluidity as in Fig. 1 of Ref. 
0). 
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Fig. 2. Quiescent thermal luminosity of several NSs in SXTs 
versus mass accretion rate compared with theoretical heating 
curves. Four ranges of (different hatching types) corre- 
spond to four scenarios of neutrino emission (Fig. 1). 




2 4 6 

lg t (yrs) 



Fig. 3. Effective surface temperatures r a °° of several isolated 
NSs versus their ages compared with theoretical cooling curves 
(from Ref. 0). Four ranges of T s °° (different hatching types) 
correspond to four scenarios of neutrino emission (Fig. 1). 

The heating curve of a low-mass NS provides an upper 
limit of whereas a heating curve of a high-mass NS 
gives a lower limit of L"^ , for a fixed scenario of neutrino 
emission. Varying the NS mass from the lower values to 
the higher we obtain a family of heating curves which fill in 
the (hatched) space between the upper and lower curves. 
In Fig. 2 we have four hatched (acceptable) regions for 
four scenarios of neutrino emission in Fig. 1. 
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4 Discussion 

These results can be compared with observations of SXTs. 
The observational data are mostly the same as in Ref. 
122- We regard as the thermal quiescent luminosity 

of NSs, and M as the time-averaged mass accretion rate. 
The value of M for KS 1731-26 is most probably an upper 
limit. No quiescent thermal emission has been detected 
|21 from SAX J1808.4-3658, and we present the upper 
limit of for this source. The limit was obtained by P. 
Stykovsky (private communication, 2004) and discussed 
in the note added in proof of Ref. [J] . The NS seems to be 
cold, but the result should be taken with caution because 
it is based on one observation (March 24, 2001) with poor 
statistics. All the data are rather uncertain and are thus 
plotted by thick dots. 

According to Fig. 2, we can treat NSs in 4U 1608-52 
and Aql X-l as low-mass NSs (the observations of Aql X-l 
can be explained [2S] assuming the presence of light ele- 
ments on the NS surface, that increases L^). NSs in Cen 
X-4 and SAX J1808. 4-3658 seem to require the enhanced 
neutrino emission and are thus more massive. The sta- 
tus of the NS in KS 1731-26 is less certain because of 
poorly determined M. Similar conclusions have been made 
with respect to some of these sources or selected groups 
(see, e.g., ^B1 _ !2E] and references therein) although four 
scenarios (i)-(4) are discussed together for the first time. 

As seen from Fig. 2, the data are still not restric- 
tive. They cannot allow us to discriminate between sce- 
narios (l)-(3) (nucleon NS cores with the direct Urea 
process, pion-condensed or kaon condensed cores). Nev- 
ertheless, the observations of Cen X-4 and especially of 
SAX J1808. 4-3658 seem to contradict scenario (4) of nu- 
cleon matter with the forbidden direct Urea process [and 
the data on SAX J1808. 4-3658 marginally contradict sce- 
nario (3)}. 

The same four scenarios can also be tested by compar- 
ing observations of thermal radiation of cooling (isolated) 
NSs with theoretical cooling curves (the dependence of the 
effective surface temperatures T s °°, redshifted for a distant 
observer, on the stellar age). These results are presented 
in Fig. 3; the data and curves are taken from Ref. [J . By 
comparing Figs. 2 and 3 we see that the data on isolated 
NSs, although more numerous, are currently less conclu- 
sive and do not discriminate between all scenarios (l)-(4) 
[although a discovery of a not too old isolated NS slightly 
colder than the Vela pulsar would contradict scenario (4)]- 

Thus, the pycnonuclear burning of highly exotic atomic 
nuclei in accreting NSs is helpful to solve the longstanding 
problem of the EOS of supranuclear matter in NS cores. 
Hopefully, new observations of SXTs and cooling isolated 
NSs will appear in the near future. Combined together, 
they will give more stringent constraints on the EOS of 
dense matter. 
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